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Abstract 

We study the processes gg — * qiq^A, with q = t,b, within the theoretical framework 
of the Supergravity (SUGRA) inspired Minimal Supersymmetric Standard Model 
(MSSM) at the leading order (LO) in perturbative Quantum Chromo-Dynamics 
(QCD). Other than constituting a novel production mechanism of the neutral, CP-odd 
Higgs particle, they also allow one to relate the size of the corresponding production 
rates to the ratio of the vacuum expectation values of the two Higgs fields tan /3, to 
the trilinear couplings Aq and to the sign of the Higgsino mass term sign(//). This 
interplay is made easier by the absence of any mixing in the q\qiA vertices, contrary 
to the case of all other Higgs bosons of the theory. 



1 Introduction and motivation 



The CP-attribute of the pseudoscalar Higgs boson induces several other oddities in its 
behaviours, with respect to all other Higgs scalars of the MSSM, that render such a particle 
a very attractive candidate for phenomenological studies. 

For example, in the qqA Feynman rule, where q represents an ordinary heavy quark 
(hereafter, q = t,b), there is no dependence on the Higgs mixing angle, a, contrary to the 
case of the CP-even scalars, H and h. As for the charged Higgs, in the qq'H^ vertex, 
another mixing, this time in the quark sector, is introduced via the Cabibbo-Kobayashi- 
Maskawa matrix. A neat consequence of this is the steep rise(fall) of the production cross 
sections of the A boson whenever this is emitted by a heavy down(up)-type quark, for 
increasing(decreasing) tan/3(cot/3) For all the other Higgs scalars, such monotonic 
behaviour is spoiled by the presence of angular terms, typically sines and cosines of a, so 
that only in extreme regions of the MSSM parameter space such peculiar dependence on f3 
can be recovered. 

If one further considers the Higgs couplings to the scalar partners of ordinary heavy 
quarks in Supersymmetric (SUSY) theories, the left- and right-handed squarks q x , with 
X = R,L, then it is easy to verify that mixing angles relating their chiral and physical mass 
eigenvalues do not enter the q~\q~iA vertex, i.e., the one involving the observable squarks 
(the subscript 1(2) referring to the lightest (heaviest) of them). Indeed, this is not the case 
for the corresponding couplings of the H, h and if ± scalars. 

One can see this in the context of SUGRA models [[|, with the minimal particle content 
typical of the MSSM (henceforth denoted as M-SUGRA, the environment we choose for 
our analysis) f||, Q , where the relevant Feynman rules for the squark-squark-Higgs vertices 
can be written in the physical basis gi )2 as follows []: 

= c q c q'^®qLq' L s q s q'^®qRq' R c q s q'^®qLq' R s q c q' ^^q R q' L ^ 
•-Wg 2 <?2 = s q s q'^$qLq' L + c q c q'^®qRq R ~ s q c q'^$qLq R ~ c q s q'^$qRq' L -> 
^*9i<?2 = ~ c q s q'^^q.Lq' L s q c q'^$qRq' R c q c q'^®qLq' R ~ s q s q'^^qRq' L i 

^*92?i = — S q C q'^q L q' L + C q S q' ^<S>q R q' R ~ S q S q'^q L q' R + C q C q'^q R q' L - W 

Here, the symbol $ denotes cumulatively the five Higgs scalars of the MSSM, $ = H,h, A, 
and H . All the \<s>q x q' , 's appearing in eq. ([!]) can be found, e.g., in the Appendix of Ref. Q. 
These are function of the five independent parameters defining the M-SUGRA model: the 
universal scalar and gaugino masses Mq and M1/2, the universal trilinear breaking terms 
A , the ratio of the vacuum expectation values (VEVs) of the two Higgs fields tan (3 = v/v', 
and the sign of the Higgsino mass term //. 

The squarks mixing angles too, i.e., 9 q , with q = t,b, can be written in terms of the 

1 Note that these Feynman rules are valid in a general SUSY model with minimal content independent 
of the origin of the soft breaking terms. 
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above four M-SUGRA parameters, as (here, s q = sin 6 q and c q = cos 6 q ) 
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with the Z-boson mass and s 2 



w = sin 2 6*v^ the (squared) Weinberg angle, m t and 
m;, the top and bottom masses, where A t and A b are the mentioned (top and bottom) 
trilinear couplings at the EW scale, while Mq 3 , M^ 3 and Mq 3 are the running soft SUSY 
breaking squark masses of the third generation, as obtained from starting their evolution 
at a common Grand Unification Theory (GUT) scale set equal to Mq. 

Now, it should be noticed that in the case of the CP-odd Higgs boson, i.e., $ = A, if 
one reverts the chirality flow in the vertex \Aq L q R , the corresponding Feynman rule changes 
its sign 0Q, 

■^Aq L q R = — ^Aq B q L , (3) 

so that, by making use of eq. (|TJ) (where A^^ = XA&fc = M), one can conclude that 
the vertices A At ~ lf ~ 2 and A^^ are independent of the mixing angles 6 t and Of,. In fact, the 
Feynman rules for those vertices reducef]to (here, = 47ra em / '$w an d = ^1(1 — s w)) 



x 9Wmt . . , 



A ^ = -|^ (/i - Afetan/5) - (4) 



These are precisely the couplings entering the processes that we are going to discuss: 



99 QiQzA 



t,b. 



(5) 



It should by now be obvious to the reader our intent in this paper. Namely, to study the 
dependence of the production rates of the scattering processes (|J) on the low-energy SUSY 
parameters, in order to pin down their actual value at the GUT scale, thus constraining 
the SUSY scenario which lies behind the MSSM, through experimental measurements of 
physical observables. 

Needless to say, such a task is greatly facilitated in case of pseudoscalar Higgs produc- 
tion, $ = A, as we have just shown that the vertex expressions for the other cases, when 
$ = H,h and H ± , are much more involved (i.e., they contain as additional free parameters 
the Higgs and squark mixing angles), so that it is inevitably much more difficult to extract 
useful information from the corresponding production rates. However, given that the final 
signatures of all possible 'gluon gluon — > squark- squark- Higgs' processes, after the decays of 
the heavy objects, are at times very similar, one cannot subtract oneself from studying the 

2 Indeed, this is the reason why vertices of the form Aq x q x i with X = x' arc prohibited at tree-level. 

3 Hcre, we neglect the CP-violating phases, by assuming that they have small values, as preferred by 
the measurements of the Electric Dipole Moments ||. The case in which such phases are sizable will be 
addressed elsewhere || . 



2 



whole of such a phenomenology. This is beyond the scope of this short note, though, and 
we will address the problem in a forthcoming publication 0. Furthermore, in that paper, 
we will also discuss more closely another relevant aspect of Higgs production in association 
with heavy squark pairs, that is, the fact that such processes can furnish additional pro- 
duction mechanisms of Higgs bosons, to be exploited in the quest for such elusive particles, 
somewhat along the lines of Ref. |7]], where the final state with both light Higgs, h, and 
stop pairs, tit*[, was considered. Finally, notice that, given the current limits on squark 
and Higgs masses f| |[, the only collider environment able to produce a statistically signi- 
ficative number of events (||) is the LHC (i/i = 14 TeV), to which we confine our analysis. 
Incidentally, at such a machine, the contribution to squark-squark-Higgs production via 
quark- ant iquark annihilations is negligible compared to the gluon-gluon induced rates p], 
so we will not consider the former here. 

The plan of this letter is as follows. The next Section briefly outlines how we have 
performed the calculations. In Sect. |3] we present and discuss our findings. The conclusions 
are in the last Section. 



2 Calculation 

The techniques adopted to calculate our processes will be described in detail in Ref. f|, 
where also the formulae necessary for the numerical computation of the Feynman ampli- 
tudes will be given. Here, we only sketch the procedure, for completeness. 

There are 10 LO Feynman diagrams for each of the two processes (H): see Fig. 1 of 
for the relevant topologies. These have been calculated analytically and integrated 
numerically over a three-body phase space. While doing so, they have been convoluted 
with gluon Parton Distribution Functions (PDFs), as provided by the LO set CTEQ(4L) 

00- 

The centre-of-mass (CM) energy at the partonic level was the scale used to evaluate 
both the PDFs and the strong coupling constant, a s . We have used the two- loop scaling for 
the latter, with all relevant thresholds [|T!| onset within the MSSM (as these are spanned 
through the Q 2 evolution of the structure functions), in order to match the procedure we 
have adopted in generating the other couplings, these also produced via the two-loop RGEs 

0- 

Depending on the relative value of the final state masses in (|5]), whether m, 2 is larger or 
smaller than (m^ + m^), the production of the pseudoscalar Higgs boson can be regarded 
as taking place either via a decay or a bremsstrahlung channel. We have treated the two 
processes on the same footing, without making any attempt to separate them, as for the 
time being we are only interested in the total production rates of the 2 — > 3 processes 
(H). In this respect, it should be mentioned that the partial widths entering in our MEs 
are significantly smaller than the total decay widths, so that processes ([5]) do retain the 

4 The systematic error due the gluon behaviour inside the proton has been investigated by comparing 
the CTEQ rates with those obtained by using other LO fits, such as MRS-LO(09A,10A,01A,07A) Q. 
Typical differences were found to be less than 15-20% §. 
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dynamics of the squark-squark-Higgs production vertices also at decay levelQ. 

Regarding the numerical values of the M-SUGRA parameters adopted in this paper, we 
have proceeded as follows. For a start, we have set M = 200 GeV and Mu 2 = 100 GeV. 
For such a choice, the M-SUGRA model predicts squark and Higgs masses in the region 
of 100-400 GeV, so that the latter can in principle materialise at LHC energies. Then, we 
have varied the trilinear soft SUSY breaking parameter Aq in a large region, (—500, 500) 
GeV, while we have spanned the tan/3 value between 2 and 45. As for /i, whereas in our 
model its magnitude is constrained, its sign is not. Thus, in all generality, we have explored 
both the possibilities sign(/x) = ±. Finally, we have gone back to consider M and M 1 / 2 in 
other mass regions. 

Starting from the five M-SUGRA parameters M , Ml/2, A , tan (3 and sign(/x), we have 
generated the spectrum of masses, widths, couplings and mixings relative to squarks and 
Higgs particles entering reactions @ by running the ISASUGRA/ISASUSY programs for M- 



SUGRA contained in the latest release of the package ISAJET [13]. The default value of 
the top mass we have used was 175 GeV. Finally, note that also typical EW parameters, 
such and sin 9w, are taken from this program. 



3 Results 

Although we will in this letter mainly concentrate on the case of CP-odd Higgs production, 
we nonetheless ought to display some typical cross sections for all processes of the form |J 

9 9 -»• Q x C* $ ' ( 6 ) 

where = t,b, = 1,2 and $ = if, h, A, H ± . This is done in Tab. I, where, for 
reference, the trilinear coupling A has been set to zero and two extremes values of tan f3, 
i.e., 2 and 40, have been selected. The corresponding mass spectrum for the particles in 
the final state of processes (|6]) is given in Tab. II. There, it is well worth noticing that 
modifying the value of tan (3 corresponds to induce quite different mass values for both 
Higgs bosons and squarks. 

From Tab. IQ, one can notice that our two processes could well yield detectable rates 
in the large tan (3 region. For a LHC running at high luminosity, some seven thousand 
such events can be produced per year. For large tan (3 values, alongside pseudoscalar Higgs 
boson production, there are at least three other mechanisms @ with observable rates, as 
one finds that, typically: 

<r(j99 tJ* 2 H; gg -> tj^h; gg -> t x t\h\ gg -> tJ* 2 A; gg -> bJ>%A) > 10~ 2 pb. (7) 

Among these, it is gg — ► tit^H that shows a strong dependence on sign(/x), whereas all 
other reactions in (0) are rather stable against variations of the latter. 

5 Generally, the dominant decay channels of squarks are to the lightest neutralino/chargino and the 
'parent' quarks (see Tab. V later on). 

6 Note that, here and in the following, the production rates do not include charge conjugation. 
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Table I: Total cross sections for processes of the type g g — > g x g'* $, where = 6, t, 
= 1,2 and $ = if, h, A, H ± 1 in the MSSM, at the leading order in perturbative QCD, 
for selected values of tan/3 and sign(/x). The other three independent parameters of the 
model have been set as: M = 200 GeV, M 1/2 = 100 GeV and A = 0. 



At low tan /3's, the cross sections for pseudoscalar Higgs boson production are presum- 
ably too poor to be of great experimental help, in both cases of sbottom and stop squark 
production. Even assuming 100 fb _1 of accumulated luminosity per year at the LHC, 
only a handful of events of the form (|5|) can be produced, if Aq pa 0, mainly through the 
gg — > t\t\A channel. (Prospects are somewhat more optimistic if the common trilinear 
coupling is much smaller than zero, but in the case of stop squarks only: we will come back 
to this point later on.) Anyhow, for small tan/? values, the three channels gg — > t x t*,h can 
boast very large production rates. Moreover, the dominant one, when x — x' — 1> exhibits 
a strong sensitivity on the sign of /i, as the two cross sections obtained for sign(/i) = ± 
differ by about an order of magnitude. 
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Therefore, both in the low and high tan (3 regime, pseudoscalar Higgs boson production 
is in general less effective than other channels in constraining the sign of the Higgsino mass 
term. This aspect will however not be investigated any further here, as it will be addressed 
in detail in forthcoming Ref. |J. In fact, here we are more concerned with the fact that 
reactions ([5]) are very sensitive to tan/5, much more than any other squark-squark-Higgs 
production channel: compare the first two columns in Tab. I with the last two, particularly 
for gg — > bxb^A (differences are over four/five orders of magnitude !). Other competitive 
mechanisms in this respect are (in the observable region, say, with a cross section above 
10~ 3 pb): gg — > t\t* 2 H , gg — > b\b\H and gg — > tib^H*. They are however suppressed, in 
general, as compared to the two CP-odd production channels, so for the time being we 
leave them aside for future studies || and concentrate exclusively on processes (j5|). 
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Table II: Masses of stop and sbottom squarks and Higgs bosons in the MSSM in the small 
and large tan /3 region, for both positive and negative values of \i and universal boundary 
conditions: M = 200 GeV, M 1/2 =100 GeV and A = 0. 



* This value is actually excluded from the latest LEP bounds on the light and CP-odd Higgs boson 
masses |J: mh,rriA ^ 80 GeV. Nonetheless we keep it here for illustrative purposes, as representative of 
the condition tua S> mh, typical of the small tan/3 regime. 



Figs. 1 and 2 further enlighten the tan (3 dependence of pseudoscalar Higgs boson pro- 
duction, in association with stop and sbottom squarks, respectively, as we have now treated 
A as a variable parameter. Indeed, the typical behaviour seen in Tab. I for A = in 
gg — > tit^A and gg — > bib^A persists for all other values of A considered. 

The variation with tan/3, and particularly the steep rise at high values of the latter, 
can be understood in the following terms. For large tan /3, the squark-squark Higgs boson 
couplings of eq. (f§) can be rewritten in the form 

A ^"~2M^' ^am 2 — — 2M^7 ' (8) 
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That is, the coupling which is associated with the sbottom pair is proportional to tan/?, 
so that, eventually, the total bib^A cross section will grow with tan 2 /? while the coupling 
related to the stop pair takes on constant values. In the latter, the enhancement of the 
tit^A cross section with tan/3 is rather a phase space effect since, as tan/? increases, the 
CP-odd Higgs boson mass decreases considerably (the squark masses changing much less 
instead), as we can see from Tab. II. Of course the same is valid in the former case and that 
is why our Figs. follow the pattern a(gg — > bib^A) > cr(gg — > tit^A) at large tan/?. 

Despite of the abundance of tit^A and bib^A events at large tan /?, overwhelming contri- 
butions involving the light stop t\ and light Higgs scalar h, i.e., t\b\h and, particularly, tit^h 
(see Tab. I), would however dominate the squark-squark-Higgs production phenomenology. 
Therefore, it might seem at first glance that reactions (|5|) cannot possibly be disentangled, 
further considering that at large tan /? the dominant decay modes of both h and A scalars 
are into bb pairs [JH| . This need not to be true though. In fact, the reader should recall two 
important aspects. Firstly, the lighter scalar quark mass, , will most likely be known 
well before a statistically significative sample of events of the form (P) can be collected. 
It follows that its knowledge can be exploited to remove gg — > q x q*/§ candidates with 
exactly two light stop squarks, thus also unwanted iit\h final states. Secondly, in the bb 
channel, one should expect experimental mass resolutions to be smaller than the typical 
mass differences tua — rrih ^ 20 GeV seen for tan (3 < 35 (see Tab. II) 0Q Needless to say, 
the light scalar h ought to have been discovered (and rrih measured) by then, for the sake of 
the all SUSY theory, so that a suitable selection of bb pairs far away from the resonance 
(or, at worse, an event counting operation in the case of overlapping h and A mass peaks, 
at extremely large tan/?) would aid to reduce t\t* 2 h events also. 

As for the low tan/3 region, as intimated a few paragraphs above, we can appreciate 
in Fig. |l] the beneficial effect of a large and negative value of Aq, in terms of the sign(/x) 
dependence of the A production rates. For example, for tan/? ~ 2, the cross sections for 
\l > and /i < take very different values, by an order of magnitude. In fact, one has 
that, for A < - 300 GeV, 

o{gg -> ht'^A) ~ 10~ 3 pb for /i > 0, a(gg -> tit* 2 A) ~ 10~ 4 pb for /i < . (9) 

Therefore, in this scenario one might aim to constrain the actual value of sign(/x) from 
the tit^A final states alone, given such a large difference. Unfortunately, the total number 
of events (after having considered the decay rates, the finite efficiency and resolution of 
experimental analyses, etc.) is again not so large, so that one would presumably be better 
off by relying on reaction gg — ► tit^h. In this respect though, one thing is worth spotting, 
i.e., the much larger value of as compared to if tan/? is small, see Tab. II. A 
consequence of this is that the decay patterns of the two Higgs bosons are very different. 
Whereas the light one would only decay into bb pairs, the pseudoscalar one would mainly 



yield tt pairs []14| . Given the huge QCD noise of the LHC, the latter might in the end 
become a competitive approach, especially if a clean electron/muon tag can be achieved in 
the (anti)top decays. 

7 The intrinsic width T$ of the h and A Higgs bosons is in that mass range about 2.5 GeV at the most. 
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But, let us now turn our attention to the other strong dependence of the production 
rates of gg — ► t\t\A and gg — > bib^A: the one on the common trilinear coupling A . This 
is in fact the most noticeable feature of both Figs. [l|-|2]: that the sensitivity to A of the 
production cross sections provides the unique possibility of constraining, possibly the sign, 
and hopefully the magnitude, of this fundamental M-SUGRA parameter. Indeed, we have 
come to believe that this is the main novelty that should be attributed to the phenomeno- 
logical potential of the processes we are studying, as one might quite rightly expect that the 
determination of tan/3 will come first from studies in the pure Higgs sector (i.e., via SM- 
like Higgs production and decay mechanisms), especially considering the theoretical upper 
limit on nih- Should this be the case, far from overshadowing the usefulness of reactions 
(§), the knowledge of tan/? would further help to constrain A Q . Let us see how. 

For a start, to observe by the thousand events involving A-production with pairs of stop 
quarks would induce the following reasoning: 

f -500 < A < for 22 < tan/5 < 40 and \l = ±. 
a(gg -> i x t* 2 A) > 10~ 2 pb { ' (10) 

{ < A < 500 for 40 < tan (3 < 42 and fi = ±. 

That is to say, unless 40 < tan (3 < 42 (a very small corner of the M-SUGRA parameter 
space), to observe such a rate of i\t\A events would mean that A Q is necessarily negative 
(whichever the sign of /i). 

Incidentally, we would like the reader to spot in Fig. |] that are the lines corresponding 
to A = 300 GeV (denoted by the arrow) those stretching to the far right of the plot, 
thus inverting the trend of decreasing rates with growing A . In other terms, such curves 
represent a true lower limit on the value of this cross section (practically for all tan /3 
values), so that the latter is bound to be in the range 

10~ 4 pb < a(gg -> U%A) < lO^ 1 pb , (11) 

values well within the reach of the LHC luminosity ! 

Similarly, one can proceed to analyse gg — ► bib^A from Fig. ^| Schematically, 

f -500 < A < for 28 < tan/3 < 40 and p = ±. 
a(gg -> hb* 2 A) > 10~ 2 pb \ (12) 

{ < A < 500 for 40 < tan/3 < 42 and \i = ±. 

Once again, to observe bib^A signals at such a rate would force A to be negative over most 
of the M-SUGRA parameter space. 

As for peculiar trends in Fig. two behaviours worth commenting on are the following. 
Firstly, that the production rates decrease with diminishing tan/3 much more than they 
do in case of stop production, particularly if tan/3 < 30. Secondly, that the cross sections 
exactly vanish in the case A = 500 GeV, when tan/3 = 24(27) if \x > 0(/i < 0), as induced 
by the oc (// — A b tan/3) dependence of the production vertex, when \p\ |/i fe tan/3| and 
At changes its sign. 

Another aspect made clear by both these two figures is that current experimental bounds 
tend to exclude only extreme parameter regions, i.e., where Aq is strongly negative and/or 
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where tan/3 is extremely high. On the one hand, LEP2 has almost exhausted its SUSY 
discovery potential, as most of the data have already been collected and/or analysed, 
whereas at Tevatron, the present ~ 120 GeV limit on the lighter stop mass is unlikely 
to be increased by the new runs to the typical values of Tab. II. On the other hand, 
the bulk of the (A , tan/3, sign(/i)) parameter space investigated here, where processes (|5|) 
could well be detected and studied at the LHC, appears in Figs. [I|-|2| far beyond the reach 
of the present colliders. Therefore, in the very short term, one should not expect that 
new experimental limits can modify drastically the look of our plots. In particular, notice 
that the presence of ti and 62 squarks in the final state of processes (||) implies that the 
corresponding production rates at Tevatron are negligible, even for optimistic luminosities, 
because of the enormous phase space suppression (see Tab. 11)0. Therefore, we believe that, 
when the LHC will start running, most of the M-SUGRA parameter space discussed here 
will still be unexplored. 

Bringing together the various results obtained so far on A , tan/3 and sign(/z), we 
attempt to summarise our findings in Tab. III. There, we list the restrictions that can in 
principle be deduced on the above three parameters by studying the two processes (H), 
assuming that none of these quantities is known beforehand. Indeed, an enormous area 
of the M-SUGRA space can be put under scrutiny, particularly involving A and tan j3. 
The prospect of the latter quantity being already known by the time gg —> tit^A and 
gg — > bib^A studies begin would be even more exciting. In such a case, a vertical line could 
be drawn in Figs. [l|-|2], so that an accurate measurement of the production cross sections 
of processes (|) would precisely pin-point the actual value of A . 

Before closing, we study the dependence of pseudoscalar Higgs boson production in asso- 
ciation with stop and sbottom squarks on the last two M-SUGRA independent parameters: 
M and Mi/2- The main effect of changing the latter is onto the masses of the final state 
scalars, through the phase space volume as well as via propagator effects in the scattering 
amplitudes. In other terms, to increase one or the other depletes the typical cross sections 
of (H), simply because the values of all m^ x and m$ get larger. Tab. IV samples such a trend 
on four among the dominant production channels, including our two reactions (|5|). As an 
example, notice that, for Mq = 300 GeV and M1/2 = 250 GeV, all the squark masses are of 
the order > 460 GeV, whereas for the heavy Higgs bosons one has that typical values are 
> 290 GeV. Not surprisingly then, among the processes in Tab. IV, for such high M and 
Mi/2 values, the only ones to survive are those involving both the lightest squark (i.e., ti) 
and the h scalar (for which one necessarily has that ^$130 GeV) [§,0]. In comparison, 
processes ([5]) are generally suppressed, as one heavy mass m 9 - 2 is always present in the final 
states and since ^ ra^. Therefore, this last exercise shows that only light M and M ly / 2 
masses (say, below 200 and 150 GeV, respectively) would possibly allow for pseudoscalar 
production to be detectable at the LHC. 

8 The gluon luminosity is much poorer too, as compared to the CERN hadron collider. 
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± 
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± 


< 1(T 2 
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40^42 


± 


< 1(T 2 
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< 40 


± 


~ 1(T 3 
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> -500 


2^-3 


+ 
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Table III: Possible restrictions on three M-SUGRA parameters derivable from studies of 
CP-odd Higgs boson production in association with stop and sbottom squarks. 



M (GeV) 


M 1/2 (GeV) 


a(99 -> iit\h) 




a{gg -> t&A) 


a(gg -> &x&*A) 


200 
200 
200 
300 


125 
150 
200 
250 


6.8 x IO" 3 
3.4 x 10~ 3 
9.7 x 10" 4 
2.6 x 10~ 4 


6.9 x IO" 2 

4.7 x 10~ 2 

1.8 x IO" 2 

5.9 x IO" 3 


4.0 x IO" 3 
1.4 x 10~ 3 

3.1 x 10" 4 
4.8 x 10~ 5 


5.2 x IO" 3 
2.4 x 10~ 3 
6.0 x 10" 4 
9.8 x 10~ 5 



Table IV: The variation of the most significant cross sections (in picobarns) of processes 
(|>D with M and Mu2- For reference, the other three M-SUGRA parameters are fixed as 

follows: A = 0, tan/? = 40 and sign(/i) = — . 



4 Summary and conclusions 

In summary, we have studied pseudoscalar Higgs boson production in association with 
stop and sbottom squarks at the LHC, in the context of the SUGRA inspired MSSM. 
Our interest in such reactions was driven by the fact that the squark-squark-Higgs vertices 
involved, other than carrying a strong dependence on three free inputs of such a model, 
i.e., A , tan/5 and sign(/x), are not affected by the presence of additional unconstrained 
parameters describing the mixing between physical and chiral squark eigenstates. 

We have found that the cross sections of such processes might be detectable both at 
low and high collider luminosity for not too small values of tan f3. Indeed, their production 
rates are strongly sensitive to the ratio the VEVs of the Higgs fields, thus possibly allowing 
one to put potent constraints on such a crucial parameter of the MSSM Higgs sector. 
Furthermore, also the trilinear coupling Aq intervenes in these events, in such a way that 
visible rates could mainly be possible if this other fundamental M-SUGRA input is negative. 
(Indeed, to know the actual value of tan j3 from other sources would further help to assess 
the magnitude of A .) As for the sign of the Higgsino mass term, sign(/i), it only marginally 
affects the phenomenology of such events. Finally, concerning the remaining two parameters 



10 



(apart from mixing effects) of the M-SUGRA scenario, i.e., Mo and M1/2, it must be said 
that their values should be such that they guarantee a rather light squark and Higgs mass 
spectrum, in order the latter to be within the reach of the LHC. 

In conclusion, we believe these processes to be potentially very helpful in putting strin- 
gent limits on several M-SUGRA parameters and we thus recommend that their subsequent 
decay and hadronic dynamics is further investigated in the context of dedicated experimen- 
tal simulations, which were clearly beyond the scope of this short letter. As a matter of 
fact, of all possible (eighteen in total) squark-squark-Higgs production modes, involving 
sbottoms, stops and all Higgs mass eigenstates, we have verified that those including the 
pseudoscalar particle are always among the dominant ones, so that one should not expect 
the presence of the others to dash away the hope of detecting and investigating the for- 
mer. In this respect, the most competing ones are those involving the lightest of the Higgs 
scalars. 



Particle BR Decay 
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A 



96% 



44% 
30% 
16% 
5% 



51% 
41% 



41% 
32% 
17% 



79% 
12% 



93% 
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92% 
6% 
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xfb 



xt b 

xt h 
biW+ 
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bb 
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T + T - 
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T + T - 
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Particle BR Decay 



Xi 



X2 



xl 



xl 



x\ 



32% 
3_2% 
14% 
10% 
10% 



Xiud) 

x\cs ] 
x o r+z/t 

XiH + v ] 



47% 

24% 
10% 
9% 



xtz 

Xth 
xtA 



69% 

10% 



x°M 



3_5% 
3_5% 
15% 



xtw- 

XiW + 

x\z 



39% 
39% 
6% 
6% 



xtw~ 

XiW + 
X\h 

x\z 



Table V: Dominant decay channels and branching ratios (BRs) of final state (s)particles 
in (g), for M = 200 GeV, M 1/2 = 100 GeV, A = 0, tan (3 = 40 and sign(/i) > 

t Via off-shell W+. * Via off-shell h. 
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This particle has however a rather different decay phenomenology from that of the CP-odd 
Higgs in most cases, whereas whenever this is not true, previous knowledge of (stop and 
sbottom) squark and/or Higgs mass values can be of some help, so that in the end it should 
not be difficult to disentangle the two scalars. 

For example, let us consider the signal gg — > tit 2 A and some possible signatures of it. 
For the choice of parameters given in the caption and in the fourth column of Tab. I, it 
yields some 3,000 events per year at the LHC. From Tab. V, one deduces that a possible 
decay chain could be the following: 



in which qq' = ud, cs and i — e, /i, r. Considering also the charge conjugated XiXi decays, 
the final signature would then be '2 jets + 46 + £^ + E m i SS \ further recalling that the 
two Xi' s an d the neutrino v produce missing energy, -E m i SS - The total BR of such decay 
sequence is 0.12 only, so that about 360 squark-squark-Higgs events would survive. One 
may further assume a reduction factor of about 0.25 because of the overall efficiency e\ to 
tag four displaced vertices (assuming e h « 0.7). This ultimately yields something less than 
100 events per year, a respectable number indeed. In addition, one should expect most of 
the signal events to lie in the detector acceptance region, since leptons and jets originate 
from decays of heavy objects. 

Such a signature has peculiar features that should help in its selection: a not too large 
hadronic multiplicity, six jets in total, each rather energetic (in fact, note that m x ± —m x o « 
30 GeV and m^m^ ^> m x ± « 63 GeV), so that their reconstruction from the detected 
tracks should be reasonably accurate; a high transverse momentum and isolated lepton 
to be used as trigger; large -E m i SS to reduce non SUSY processes; four tagged 6-jets that 
can be exploited to suppress the i W ± + light jet' background from QCD with one bb 
pair resonating at the A mass (which is well above the h one: see Tab. II and recall the 
discussion in Sect. [3] about the interplay between gg — > tit^A and gg — > tit\h events). Even 
the background from gg — > t\t\Z events, with Z — > bb, potentially very dangerous because 
irreducible and since m t ~ 2 — mi ^> M z (Tab. II), should easily be dealt with. In fact, notice 
that rriA 3> M z , so that to select only events for which M bb 96 M z would presumably allow 
one to reduce also such a noise to manageable levels. 

Notice that the M-SUGRA point just discussed corresponds to a rather low lightest 
chargino mass, though still roughly consistent with the latest bounds drawn by the Par- 
ticle Data Group (PDG) Hl5fl . However, preliminary results from LEP and Tevatron have 
meanwhile increased the limit on m+ , up to 80-90 GeV or so. Thus, we also have consid- 
ered a second parameter combination yielding sizable production rates, but now satisfying 
the latter constraint: e.g., that in the first line of Tab. IV (see the caption for the parame- 
ter setup), for which m x ± = 89 GeV, right at the edge of the exclusion band. Considering 
again '2 jets + 46 + £ ± + E miss ' decays, starting from some 400 signal events every 100 




Xi + b Xi + b 

i i 

q + q' + Xi fT + v + X °i 



A 

i 

b + b 
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inverse femtobarns produced in the 2 — > 3 scattering (fifth column in Tab. IV), one ends up 
with 12 events, as the total decay BR is basically the same as before and including again 
the factor e\. The number is reduced by an order of magnitude, but still sizable. 

As a matter of fact, other signatures can possibly be even more accessible. Let us 
now take, e.g., M = M 1/2 = 125(130) [140] GeV (with again tan/? = 40, A = and 
sign(/i) = — ). For such settings, the lightest chargino mass is m x ± = 90(93) [101] GeV. 

In correspondence, one gets a{gg -> tJ^A) = 96(79) [39] fb, i.e., some 9,600(7,900) [3,900] 
events per luminosity year. For these last three combinations of M-SUGRA parameters, it 
turns out that an interesting decay sequence could be the following: 



Apart from the BR of the channel i% — > Xib, which ranges at 25%, the others are all 
dominant and close to unity. Here, the final signature is '46 + r + r _ + i? m i SS ' with a total 
BR of about 23% in all cases. Thus, after multiplying by ef, one finally gets 528(454) [218] 
detectable events every 100 fb -1 . 

This additional channel appears particularly neat thanks to an even smaller jet mul- 
tiplicity. In addition, all such jets should be rather energetic, as m x ± — m f rs 19(17) [15] 
GeV and m? — m x o « 21(25) [31] GeV. Standard backgrounds from L Z + jet' production 
could be strongly suppressed because of the absence of light-quark jets and the presence of 
four heavy ones. From one pair of these, one could further attempt to reconstruct the A 
mass, at around 114(120) [135] GeV. Finally, the large amount of E m[ss building up because 
of the four neutrinos and two neutralinos could prove to be a further good handle against 
non-SUSY processes. As for irreducible SUSY backgrounds, notice the poor decay rate 
BR(t* -> t\Z) « 7% ! (Typical stop masses are around 380(388) [406] and 240(248) [265] 
and for i 2 and t\, respectively.) 

These are just a few illustrative examples of some possible manifestations of squark- 
squark-Higgs events at the LHC. Dedicated analyses of all mechanisms of the form gg — > 
q x q*,&, for any possible combination of q(> — t, b, x^ — 1> 2 and $ = H,h, A, if ± , of their 
interplay and a simulation of possible backgrounds and detection strategies is now under 
way §. 
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Figure 1: Total cross sections for process gg — > t\t\A as a function of tan/3. These are 
plotted for both positive (dotdash) and negative (solid) values of /i and for a selection 
of A values. (The other M-SUGRA parameters have been chosen as M = 200 GeV 
and Mi/2 — 100 GeV.) The symbol " o (x) " is used to indicate data points forbidden 
by experimental bounds from direct searches of the lightest Higgs(squark) scalar, yielding 
2 < tan/? < 40(771^ > 120 GeV). Two points with the masses of the lightest stop and the 
CP-odd Higgs are also shown. 
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Figure 2: Same as in Fig. [T] for process gg — > bib^A. 
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